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Abstract. The origin of subluminous B stars is still an unsolved problem in stellar evolution. Single star as 
well as close binary evolution scenarios have been invoked but until now have met with little success. We have 
carried out a small survey of spectroscopic binary candidates (19 systems consisting of an sdB star and late type 
companion) with the Planetary Camera of the WFPC2 onboard Hubble Space Telescope to test these scenarios. 
Monte Carlo simulations indicate that by imaging the programme stars in the R-band about one third of the 
sample (6-7 stars) should be resolved at a limiting angular resolution of O'.'l if they have linear separations like 
main sequence stars ("single star evolution"). None should be resolvable if all systems were produced by close 
binary evolution. In addition we expect three triple systems to be present in our sample. Most of these, if not all, 
should be resolvable. Components were resolved in 6 systems with separations between 0'.'2 and 4'.'5. However, 
only in the two systems TON 139 and PC 1718+519 (separations 0'.'32 and 0'.'24, respectively) do the magnitudes 
of the resolved components match the expectations from the deconvolution of the spectral energy distribution. 
These two stars could be physical binaries whereas in the other cases the nearby star may be a chance projection 
or a third component. Radial velocity measurements indicate that the resolved system TON 139 is a triple system, 
with the sdB having a close companion that does not contribute detectably to the integrated light of the system. 
Radial velocity information for the second resolved system, PC 1718+519, is insufficient. Assuming that it is not 
a triple system, it would be the only resolved system in our sample. Accordingly the success rate would be only 
5% which is clearly below the prediction for single star evolution. We conclude that the distribution of separations 
of sdB binaries deviates strongly from that of normal stars. Our results add further evidence that close binary 
evolution is fundamental for the evolution of sdB stars. 
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1. Introduction 

Subluminous B (sdB) stars dominate the populations of 
faint blue stars of our own Galaxy and are found in 
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both the disk (field sdBs) and globular clusters (Moehler 
et al. 1997). Observations of elliptical galaxies with the 
Ultraviolet Imaging Telescope (Brown et al. 1997| ) and the 
Hubble Space Telescope (Brown et al. 2000) have shown 
that these stars are sufficiently common to be the dom- 
inant source for the "UV upturn phenomenon" observed 
in elliptical galaxies and galaxy bulges (see also Greggio 
1999| ). Their space distribution and kine- 



Renzini 1990 



matical properties indicate that the field stars belong to 
the intermediate to old disk population (de Boer et al. 



1997; Altmann & de Boer 200C) 



However, important questions remain concerning their 
formation process and the appropriate evolutionary 
timescales. This is a major drawback for the calibration 
of the observed ultraviolet upturn in elliptical galaxies as 
an age indicator. 

It is now generally accepted that the sdB stars can 
be identified with models for Extreme Horizontal Branch 
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(EHB) stars burning He in their core, but with a very 
tiny (<2% by mass) inert hydrogen envelope (Heber 1986; 
Saffer et al. 1994). An EHB star bears great resemblance 



to a helium main-sequence star of half a solar mass and its 
further evolution should proceed similarly (i.e. directly to 
the white dwarf graveyard) as confirmed by evolutionary 
calculations (Dorman et al. 1993). 

How stars evolve to the EHB configuration is contro- 
versial. The problem is how the mass loss mechanism in 
the progenitor manages to remove all but a tiny fraction 
of the hydrogen envelope at precisely the same time as 
the He core has attained the minimum mass (« 0.5Mq) 
required for the He flash. 

Both non-interacting (scenario i) , and interacting (sce- 
narios ii and iii) evolutionary scenarios have been proposed 
to ex plain the origin of the sdB stars (see Bailyn et al. 



1992) 



(i) Enhanced mass loss on the red giant branch (RGB) be- 
fore or during the core helium flash may remove almost the 
entire hydrogen-rich envelope. This is usually modelled by 
increasing the 77 factor in the Reimers (1975) formula to 
estimate mass loss rates for RGB stars. It has been con- 
jectured that the mass loss rates increase with increasing 
metallicity, implying that metal rich populations should 
produce more sdB stars than metal poor ones. Birthrate 
estimates for sdB stars indicate that only 2% (Heber 1986) 
or even less (0.25% to 1%, Saffer & Liebert |1995| ) of the 
RGB stars need to experience such enhanced mass loss. 
Evidence that this is possible comes from the existence 
of RR Lyrae stars of population I which must also have 
lost half of their mass during evolution. In both cases the 
physical reason for such strong mass loss is not yet under- 
stood. 

(ii) Mengel et al. ( 1976| ) suggest that sdBs could be formed 
from binaries in which mass transfer starts on the red 
giant branch and results in a reduction of the hydrogen 
envelope prior to the helium core flash. Hence all sdBs 
star are predicted to be found in close binary systems. 

(iii) An alternative scenario was proposed by Iben ( 1990| ), 
who pointed out that sdBs can be formed from mergers 
of helium white dwarf binary systems. Iben & Tutukov 
( |l992| ) estimate that 80% of the sdBs could have been 
formed by mergers. Hence the frequency of sdBs still being 
in binaries should be at most 20%. 

Several dozens of objects with composite spectra con- 
sisting of an sdB and a dwarf G-K star have been discov- 
ered (e.g. Fe rguso n et al. |1984 Theissen et al. |1993| , |1995 ; 
Allard et al. 1994) which implies that the binary frequency 
of sdBs is 50% or more (Allard et al. 1994). The observed 
large binary frequency rules out the merger scenario (iii) 
and we are left with scenarios (i) and (ii), i.e. either the 
sdB binaries are mostly wide systems that did not inter- 
act so that the sdB precursors have evolved independently 
from the companion (i), or they are close systems formed 
by interaction of the sdB precursor with the companion 
star (mass exchange, ii). 

The high spatial resolution of the Planetary Camera 
(PC) on board the Hubble Space Telescope (HST) allows 



to perform a crucial test. As we will show in this paper, 
it should be possible to resolve a significant fraction of 
the known composite spectrum systems containing an sdB 
star if scenario (i) is correct, i.e. if the systems have a dis- 
tribution of separations like normal main sequence bina- 
ries (Duquennoy & Mayor 1991). The interacting scenario 
(ii), however, predicts that all sdB stars reside in short 
period (P<100d) binaries and consequently none of the 
systems should be resolvable even with the PC. In order 
to measure their distribution of separations we have im- 
aged 23 sdB binary candidates with the PC by taking 
advantage of the snap shot mode of HST observations. 

2. Observations and Data Analysis 

2.1. Target selection and optical spectroscopy 

For the snapshot observations a target list of fifty of the 
brightest sdB star binary candidates was extracted from 
an updated version of the Kilkenny et al. (198S) cata- 
logue, supplemented by two stars which we discovered 
in the course of follow-up spectroscopy of hot stars from 



the Hamburg-ESO survey (see Edelmann et al. 2001a). 23 
stars from this target list were actually observed with the 
Wide Field Planetary Camera 2 (WFPC2) onboard the 
HST during our snapshot project, i.e. they were sched- 
uled for observation to fill small gaps in the HST sched- 
ule. All stars have published photometry (see Tables B.l 



and B^), but only 16 have published optical spectroscopy. 
Therefore additional spectra were obtained at the Calar 
Alto and ESO observatories (see App endix ^ f or de tails 
and plots of the spect ra in Fig. AT and Fig. |A.2| ). As 
can be seen from Fig. |A.l| spectral features (Cai, Can, 
Mgi and/or Fei) indicative of a cool star are clearly 
present in the spectra of PC 1309-078, PC 0942-^461, 
HE 0430-2457, HE 2213-2212, and PG 2148-^095 in ad- 
dition to the Balmer and helium lines of the sdB. Hence 
these objects are spectroscopic binaries consisting of an 
sdB star and a cool companion. PG 0942-1-461 has already 
been observed by Mitchell (1998), who, however, did not 
note the binary nature of the star. We do not find any evi- 
dence for a cool companion in the spectra of the sdB stars 
PG 1558-087 and KPD 2215-^5037 (see Fig. ^). We 
also re-analysed a published spectrum of PG 2259-1-134 
(Theissen et al., 1993) and do not find any spectroscopic 



evidence for a cool companion. PG 0105-1-276 turns out to 
be not an sdB star but a helium-rich sdO star and does 
not show any spectroscopic evidence for a cool companion. 
Therefore our sample consists of 19 composite spectrum 
objects plus four stars showing only photometric evidence 
for a companion. One of these four stars (PG 0105-1-276) 
also does not belong to the programme sample since it is 
an sdO star. 

2.2. WFPC2 data 

We observed the candidate binary systems with the PC 
chip of the WFPC2. If the cool companion is a main 
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Table 1. Programme stars: coordinates, observation dates, and exposure times for the WFPC2 and references for the 
spectroscopic classification observations 



star 



Q1950 



1 



obs. exp. reference 













date 


time 

[s] 


PB 6107 


00''39'"31' 


+04°53'17" 


118-59 


-57?64 


990627 


3.5 


PHL 1079 


0l''35"M8= 


+03°23'00" 


144? 96 


-57?22 


981204 


4 


HE 0430-2457 


04''30'''59' 


-24°57'37" 


223-49 


-40? 55 


980417 


8 


PG 0749+658 


07'M9'"39' 


+65°50'13" 


150-44 


+30? 99 


990329 


1.8 


PG 0942+461 


09''42'"02= 


+46°08'38" 


173-11 


+48? 89 


980530 


10 


TON 1281 


10''40'''57' 


+23°24'55" 


213-62 


+60? 89 


990623 


5 


TON 139 


12''53'''39' 


+28°23'31" 


77?21 


+88? 57 


980103 


1.8 


PG 1309-078 


13''09'"09' 


-07°49'18" 


311-60 


+54? 44 


980505 


8 


PG 1421+345 


14''21'"29' 


+34°27'53" 


58-36 


+69?01 


990605 


14 


PG 1449+653 


14^49^42= 


+65°17'58" 


104^84 


+47?63 


990619 


7 


PG 1511+624 


15''11'"25= 


+62°21'00" 


99? 21 


+47?96 


990513 


14 


PG 1601+145 


16''01'"47' 


+14°32'58" 


27?15 


+43? 51 


000613 


12 


PG 1636+104 


16''36'"40'= 


+10°24'54" 


27? 00 


+34?04 


000612 


8 


TON 264 


16''47'"05' 


+25°15'13" 


45? 16 


+37?12 


990529 


10 


PG 1656+213 


16''56'"12= 


+21°15'05" 


41? 25 


+33?90 


980301 


12 


PG 1718+519 


17'' 18"" 35= 


+51°55'05" 


79? 00 


+34?94 


990427 


7 


PG 2148+095 


2l''48"'41= 


+09°30'39" 


66? 78 


-32? 84 


990411 


4 


HE 2213-2212 


22''13"'38= 


-22°12'26" 


32? 63 


-54? 50 


981207 


8 


BD -7°5977 


23'' 15"" 12= 


-06°44'56" 


71? 55 


-59?65 


981125 


0.3 



Moehler et al. (09901) 
Theissen et al. ( [l995| ) 
this paper 



Saffer (1991) 



Heber et al. ( |1991| ) 

Jeffery fc PoUacco (|l998|) 



Green (1997) 



(199C) 



Ferguson et al. (1984) 
Ferguson et al. (1984) 
Moehler et al. ( [l990| ) 
Moehler et al. 
Ferguson et al 
Ferguson et al 
Theissen et al. 
Ferguson et al 
Theissen et al. 
this paper 

this paper 

Viton et al. ([f99l|) 



(1984) 
(1984) 
(1993) 
(1984) 
(L995) 



stars without spectroscopic evidence for a cool companion 



PG 0105+276 
PG 1558-007 
KPD 2215+5037 
PG 2259+134 



01''05"'32= 
15''58'"39= 
22''15"'25= 
22''59'"16= 



+27°36'53' 
-00°43'26' 
+50°37'48' 
+13°22'31' 



127?46 -34?84 980226 14 this paper, new type: He-sdO 

9?34 +36?51 990424 7 this paper 

99?71 -4?91 961213 7 this paper 

86?36 -41?31 000615 10 Theissen et al. ( [l993| ), this paper 



sequence star, both components should be of compara- 
ble brightness in the R band and we therefore used the 
F675W filter of the WFPC2. We obtained four observa- 
tions of each target, which were offset relative to the first 
one by (-11,-5.5), (-16.5,-16.5), (-5.5,-11) pixels. We 
first rebinned the data linearly to a step size of 0.5 pix- 
els and then aligned them according to the offset pattern 
mentioned above. We then determined the median value 
of the four aligned images to avoid cosmic ray hits and 
hot pixels and used these median-averaged images for vi- 
sual inspection. All flux measurements are performed on 
manually cleaned average images to ensure proper flux 
conservation. 

The median-averaged images were first inspected by 
eye to see if any companion could be detected. Only 6 stars 
(cf. Fig. 1^) showed obvious nearby stars and the angu- 
lar separations and brightness differences can be found in 
Table ^. The brightness differences were determined using 
the command INTEGRATE/APERTURE from MIDAS, which 
performs an aperture photometry with a given radius. 
Aperture photometry is difficult for TON 1281, TON 139, 
and PG 1718+519, due to the small distance of the com- 
ponents. The sky background was determined in an empty 
region using the same aperture as for the stars. 

To get a more quantitative estimate of possible com- 
panions we fitted two-dimensional Gaussians with variable 



Table 2. Separation and estimated brightness differences 
for the components of the 6 resolved binaries. The photo- 
metric data available for HE 0430—2457 do not allow to 
estimate a temperature or distance of the sdB. 



system 


separation 


brightness 




angular 


linear 


diflterence 






[AU] 


AF675W 


PG 0105+276 


3'.' 37 


3700 


0-^9 




4'.' 48 


4900 


i-^e 


HE 0430-2457 


1'.'25 




2-^1 


TON 1281 


0'.'22 


250 


3-^7 


TON 139 


0'.'32 


300 


0-^8 


PG 1558-007 


2'.' 80 


2500 


3';'1 


PG 1718+519 


0'.'24 


230 


0-^8 



angle of the major axis to all shifted and co-added tar- 
get images and compared the results to fits obtained for 
archive point-spread functions (PSFs; F675W filter, PC 
chip). The archive PSFs define a good correlation between 
the length of the two axes, which is shared by most tar- 
get PSFs (see Fig. H). Besides the resolved binaries (where 
stray light can affect the determination of the axis ratio) 
four stars deviate from the main correlation between ma- 
jor and minor axis (see Fig. |): PG 2148+095 (2.03/1.26), 
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PG0105+276 





HE0430-2457 



Fig. 1. The images of the resolved binaries. The bar in each image corresponds to 1'' 



KPD 2215+5037 (2.38/1.61), TON 264 (2.35/1.83), and 
PG 0749+658 (2.36/1.87). 



We used DAOPHOT (Stetson |1987D to obtain an average 
PSF from those target stars that share the axis-relation 
of the archive PSFs. This "target PSF" was then used 
to deconvolve all systems that are either resolved by eye 
or show deviations from the axis-relation defined by the 
archive PSFs. No additional components were resolved in 
this process, but we could verify the brightness differences 
between the components of the resolved systems listed in 
Table ^, which were reproduced by DAOPHOT also for small 
separations. 

For 13 of our target stars a homogeneous set of gro und- 
based Re measurements exist (Allard et al. 1994, see 
Table |B.2| ). Comparing those data to the instrumental 
^67514^ magnitudes integrated within an aperture of 0'.'5 
radius 



flux J/ ^ 

i^675W^ -2.51og 



exposure time 



we find that most of the stars lie on a line with slope 1 (ex- 
cept KPD 2215+5037 and PG 1601+145, see Fig. |). From 
the 11 stars on the line we determine a zeropoint of 21 "^21 
± 0"^02. From the WFPC2 data handbook we determine 
a zeropoint of 2l"^9 (gain 14, including an aperture cor- 
rection of — 0™1) that has to be corrected to Cousins R 
by adding — 0™65 (assuming a spectral type of A5 for the 
combined spectra of our binary stars), yielding a final ze- 
ropoint of 21™25, in agreement with our empirically deter- 



'a 



s 1.5 



& ■» archive PSF 
O • target PSF 



O 



o 



4^ A'^ oo 

o 



4- ^ 



1.5 2 

major axis [PC pixel] 



2.5 



Fig. 2. The major and minor axes of the point spread 
functions for the target stars (circles, filled symbols mark 
brightest star of resolved binaries) and of archive point- 
spread functions (triangles, filled ones mark stars with po- 
sitions on the PC chip close to our targets). 



mined zeropoint. Since our empirically derived zeropoint 
automatically takes into account the unusual flux distri- 
bution of our binary stars we decided to use it to calculate 
Rest given in Table B.2. 
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PG2148+095 



KPD2215+5037 



PG1421+345 



T0N264 



PG0749+658 



PB6107 



Fig. 3. The images of the unresolved stars (PG 2148+095, KPD 2215+5037, TON 264, PG 0749+658) which show 
deviations from the standard PSF shape (see text). The images of PB 6107 and PG 1421+345 are well matched by 
the standard PSF shape and are displayed for comparison. Note that - in contrast to all other stars displayed here - 
there is no spectroscopic evidence for binarity of KPD 2215+5037. The bar in each image corresponds to 1". 



CD 



-7.5 



-8.5 




21.ai±0.02 



14.5 



Fig. 4. The instrumental F675M^ magn itudes compared 
to the Rc data from AUard et al. (1994). The open sym- 
bols are KPD 2215+5037 and PG 1601+145. The line 

flux „ -sky , 



marks the relation Rc — —2.5 log 



exposure time 



-21.21. 



3. Spectral energy distribution 

To obtain an upper limit to our resolution we tried to es- 
timate the R brightness of the cool companion by fitting 
the available photometric data of those stars that have 
sufficient measurements. In order to disentangle the flux 
of the hot star from that of the cool star we analyse the 
composite spectral energy distribution. For this purpose 
ultraviolet, optical and infrared (spectro-) photometry is 
collected from literature and archives (lUE, 2MASS). To 
determine the contribution of the hot star we fit synthetic 
spectra (Kurucz, 1992| ) to the bluest part of the observed 
spectral range, i.e. lUE data plus u or u/U plus v/ B (if no 
UV data were available) and determine the effective tem- 
perature of the sdB star. In doing so we assume that the 
companion does not contribute to the flux in this wave- 
length range (cf. Fig. ||). While this is probably true for 
the lUE data, some contamination may be present in the 
u/U- and w/i?-band and consequently the temperature 
determination for the sdB star can be compromised. 

However, for some stars photometric data are so 
incomplete that no meaningful flt can be obtained. 
Aside from the FQlbW measurements discussed here 
PG 0942+461 and HE 2213-2212 have only JHK pho- 
tometry from 2MASS, which are insufficient for a fit. 
While HE 0430-2457 has BVR photometry it is still not 
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Table 3. Flux for a star with to a = 
taken from Lamia ( 



1982 



0. The data are 
p. 82 BVRcIc), 



p. 59, uvhy 

Zombeck ( |1990| , JHKutqs) and from the 2MASS Team 
(priv.comm., JHK2mass)- 



filter 


flux 


Ac 




[erg/(cm^ s A)] 


[A] 


u 


1.169 • 10"** 


3500 


V 


8.444 • 10"^ 


4110 


b 


5.826 • 10"" 


4670 


y 


3.700 ■ 10"^ 


5470 


u 


4.187 • 10"^ 


3600 


B 


6.597 • 10"^ 


4400 


V 


3.607 • 10"^ 


5500 


Rc 


2.254 • lO"'' 


6400 


Ic 


1.196 • 10"" 


7900 


J2MASS 


2.91 ■ 10"i° 


12510 


-ff2MASS 


1.11 • 10"i° 


16280 


-?f2MASS 


3.83 • 10"" 


22030 


>7uT98 


3.18 ■ 10-1° 


12500 




1.18 ■ 10"" 


16500 




4.17 • 10^" 


22000 



possible to constrain the sdB star's temperature with these 
data as i? — F is insensitive to T^s at sdB temperatures. 
To convert the magnitudes into flux values we used the 
data given in Table |. 

By comparing the measured flux in the R band to the 
model flux of the sdB star we derive the flux ratio of the 
hot vs. the cool star in the system. For those systems which 
should have a rather bright companion according to their 
photometric data we verified the flux ratio in R between 
sdB and cool companion from two c olour diagrams simi- 
lar to those used by Ferguson et al. (1984), which is best 
suited for components of co mparable brightness (for de- 
tails see Ferguson et al. 1984 ). With this method we found 
that the companion of TON 1281 is bright enough to af- 
fect also the u filter, yielding a temperature of 25,000 K 
to 27,000 K for the sdB instead of the 22,000 K given in 
Table ^ and a brightness difference AR of 0™2 to — 0™1. 
Also for PG 1601-1-345 we find a much smaller brightness 
difference (0'"l) and higher temperature (29,500 K) from 
this method than from our photometric fits. In this case 
the B filter is already affected by the cool companion. For 
reasons of consistency we keep the values from the photo- 
metric fits for these two stars in Table |[ For all other stars 
with brightness differences <0™8 the results from both 
methods were the same. To correct for interstellar red- 
dening we used the reddening-to-infinity maps of Schlegel 
et al. (1998) which give somewhat higher values than the 
older data of Burstein & Heiles ( |1982D . KPD 2215-1-5037, 
PG 1558-007, and PG 2259+134 aU lie in regi ons of quite 
high reddening according to Schlegel et al. (1998) and 



show no spectroscopic evidence for a cool companion (see 
Appendix ^). The observed apparent infrared excess can 
be explained by high interstellar reddening alone, without 
invoking the presence of a cool companion. We also find 
no evidence for a companion from available photometry 



of PG 1656+213, although there is spectroscopic evidence 
(Ferguson et al, 1984). However there are are no flux mea- 
surements redwards of V available and B and V fluxes are 
inconsistent. Therefore we keep PG 1656+213 as a pro- 
gramme star. 

Aznar Cuadrado & Jeffery (2001) present an extensive 
discussion of sdB parameters derived from energy distri- 
butions, which also includes some of the stars discussed in 
this paper. In Table |5| we present the temperatures given 
in their paper and other values collected from literature 
in comparison to the ones derived here. As can be seen 
from Table ^ differences of ±10% in Toff between different 
authors are quite common. 

The temperatures derived from the photometric data 
and from line profile fits for the stars in regions with high 
reddening agree moderately well (compare Tables ^ and 
A.l). The discrepancies may be due to small scale vari- 
ations in reddening that affect the temperatures derived 
from photometry but not those derived from line profile 
fits. 

From the photometric fit we can derive the apparent 
R magnitudes of the sdB and of the cool star and correct 
both for interstellar extinction. The uncertainty in Toff of 
about ±10% evident from Table | causes an estimated un- 
certainty in the derived brightness for both components of 
±0"^2. Knowing the absolute R magnitude of the sdB stars 
then allows to determine their distance. We use the mean 
Mv derived by Moehler et al. (1997) for hot subdwarfs in 



the globular cluster NGC 6752. They found two groups 
of hot subdwarfs, a cooler one with a mean effective tem- 
perature of 22,000K and < My > = 3'"2 (5 stars), and a 
hotter one with <Toff> = 29,000 K and < My > = 4'°2 



(12 stars). From Kurucz (|1992|) model atmospheres for 
[M/H] = we find F - = -0"l20 for Toff = 22,000K 
and -0'^'152 for 29,000 K. We therefore use Mr ^ 3'?3 for 
stars cooler than 25,000 K and Mr — 4'"4 for hotter stars. 

Using the archive point spread functions we estimated 
the minimum separation that we can resolve for a given 
brightness difference by adding two PSFs with a defined 
brightness difference and angular separation and examin- 
ing the resulting image by eye. We find the following res- 
olution limits: Aaiim (AR) = 0'.'2 (2""0), O'.'l (l^'S), 0'.'07 
(l'°0), 0'.'05 (0™5). Using the distances determined above 
we can now derive upper limits for the linear separation of 
the unresolved binaries (cf Table ^) , ranging from 50 AU 
to 210 AU. 

Table ^ shows that the brightness differences between 
the components in TON 1281 and HE 0430-2457 are 
too large to reproduce the spectral energy distribution of 
TON 1281 and the photometry of HE 0430-2457, respec- 
tively. The large brightness difference of 3™1 (from the 
WFPC2 data) for PG 1558-007 agrees with the lack of 
photometric and spectroscopic evidence for a companion. 
In the remaining two cases (PG 1718+519, TON 139) the 
brightness differences in Table ^ are somewhat larger than 
those derived from the spectral energy distribution. To see 
whether we can in principle accommodate the HST obser- 
vations by fits to the photometric data we repeated the 
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Table 4. Estimated temperature of sdB stars, resulting reddening-free brightness of subdwarf B star (i?sdB,o) and 
companion (i?comp,o)j distance d, brightness difference AR, and upper limit for linear separation aiim derived from 
upper limit of angular separation Aaum- The reddening estimates are from the maps of Schlegel et al. (1998) and we 
used An — 2.6-Eb-v- The three different temperatures for PG 1511+624 result from the three available SWP spectra. 
If no evidence for a companion can be found from available photometry no entry is given in column 4. 
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Table 5. Effective tempera tures for sdB stars derived fr om en ergy distributions by v arious auth ors. T he sources are 
Aznar Cuadrado & Jeffery ( |200l| , ACJOl), Allard et al. ( |1994 A94), Theissen et al. (|1993| , T93; |1995| , T95), UUa & 
ThejU (|1998|, UT98). 
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fits, this time enforcing the brightness difference in the R 
band obtained from the HST data. The results are shown 
in Fig. 1^ (in comparison to the original fits). Obviously the 
companion of PG 1718+519 is sufficiently bright to affect 
also the u filter, thereby rendering our assumption that 
this filter is unaffected by the cool companion obsolete. 
The fits for TON 139 do not show much difference. We 



conclude that the spectral energy distribution of TON 139 
and PG 1718+519 are consistent with the R band flux ra- 
tio measured with the HST WFPC2 camera. 
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Fig. 5. Fits of ATLAS9 model spectra (Kurucz [l99^ , [M/H] = 0) to the photometric data of PG 1718+519 (left panel, 
including lUE spectra) and TON 139 (right panel). The upper panels show the fits obtained assuming that the bluest 
photometric data points (lUE spectra and u for PG 1718+519, u and v for TON 139) are not affected by the cool 
companion. The lower panels show fits that reproduce the brightness differences measured on the WFPC2 images. 



3.1. The sdO star PG 0105+276 

Since the He-sdO PG 0105+276 does not belong to the 
programme sample, we discuss it separately. It is the only 
programme star that is resolved into three components. 
However, the two companions are quite distant from the 
primary (3'.'37 and 4'.'48, respectively). The hght of these 
companions can explain at least qualitatively the IR ex- 
cess observed by ground based aperture photometry. The 
spectrum of PG 0105+276, however, does not show any 
signature of a cool companion, probably because due to 
the orientation and the small width of the slit no light 
of the distant companions was included. The diaphragm 
used in the photometry was large (18") and included the 
companions' light. 

The brightness differences measured on the WFPC2 
image (0'?9, 1"6) for PG 0105+276 are smaller than the 
one derived from the photometric fit (1™4), i.e. one com- 
panion is brighter than expected. However, as discussed 
in Appendix 0, the true temperature (from line profile 
fitting) is much higher than the one obtained from the 



spectral energy distribution (63,000 K vs. 35,000 K) mak- 
ing the companion's luminosity obtained from photometry 
a lower limit only. 



4. Simulation of separability in binary systems 

In order to interpret our results with respect to the dif- 
ferent evolutionary scenarios we simulate binary systems 
containing main sequence (MS) companions and sdBs with 
period distributions found for normal main sequence bina- 
ries (Duquennoy & Mayor 1991[ ). Assuming that the sdB 
mass is 0.5 Mq and the MS companion mass is 1 Mq we 
convert the period distribution published by Duquennoy 
& Mayor (1991) to physical separations using Kepler's 
Harmonic law. The orientation of the axis of the system 
is then chosen to be random in space and the projected 
separation, or asini, is calculated, given the distance to 
the system which is found from the apparent and abso- 
lute brightness of the system. The orbits are assumed to 
be circular. 
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Based on the spectroscopic distances derived above 
(see Table |) we then simulate a huge number of 
such binary systems. For three stars (HE 0430—2457, 
PG 0942+461, and HE 2213-2212) the magnitude ratio 
of the components could not be determined and therefore 
the distances are unknown. We adopted the mean value of 
the other stars (Ai? = 1™1), which is consistent with their 



spectral appearance (see Fig. A.l). The numerical simula- 
tion predicts a mean value of asini = 0'.'04 and that, out 
of the 19 observed systems, we should resolve six systems 
at a resolution limit of O'.'l, one of which should show a 
separation greater than I'.'O. 

Since the orbital motion for an eccentric orbit is lower 
during phases of large separation, the time averaged dis- 
tance is larger than the semi major axis. Thus eccen- 
tric orbits would increase the detectability. Duqucnnoy & 
Mayor (1991) also provide a distribution of ellipticities 
for normal stars. If the sdB systems did not experience 
phases of binary interaction, the distribution of eccentric- 
ity should correspond to that of normal stars. We used 
Duqucnnoy & Mayor's distribution corrected for selection 
effects. For each eccentricity the ratio of the time averaged 
distance to a was calculated and finally the mean over the 
Duqucnnoy & Mayor distribution computed. We find the 
average distance of the companions to increase by 17%. 
Another mechanism that tends to increase the separation 
of the components in a sdB binary is mass loss during post- 
main sequence evolution in order to reduce the mass of the 
sdB progenitor to its present value of half a solar mass. 
Assuming that the sdB evolved from a 1 Mq main se- 
quence progenitor it must have lost 0.5 M© due to a stellar 
wind during its post-main sequence evolution. Assuming 
that the wind emanates in a spherical symmetric manner 
and does not interact with the companion the increase in 
separati on can be calculated according to ^ — 



(Pringle 1985 ), with a being the separation and Ms and 
Mc the masses of the sdB progenitor and that of the cool 
star, respectively. As a result the separation increases by 
33%. 

We repeated the Monte Carlo simulations for increased 
separations. Even when we consider both elliptical orbits 
and evolution of the orbits due to a stellar wind as de- 
scribed above the prediction increased only slightly 7 re- 
solvable stars in our sample. 

Hence we predict that 6 to 7 stars should be resolvable 
in our sample if the systems have separations consistent 
with the Duqucnnoy & Mayor (1991) distribution. 



5. Chance projections and triple systems 

In the vicinity of five programme stars we found an addi- 
tional object within a radius of 3'.'O0. We have demon- 
strated above that only in two cases (TON 139 and 
PG 1718-1-519) the relative brightnesses are consistent 



with the expectations from the deconvolution of the spec- 
tral energy distribution. The remaining three cases must 
then be chance projections or triple systems. Since the 
programme stars lie at high galactic latitudes (except 
KPD 2215-f5037, see Table we expect chance coin- 
cidences to be rare. Indeed, we do not find any additional 
object in the PC field (40"x40") except for the low galac- 
tic latitude object KPD 2215+5037. 



According to Abt & Levy (1976) 16% of multiple sys- 
tems of normal stars are triples. If the fraction of triple 
systems is the same for our sample, we expect three pro- 
gramme stars to be triple. Most of these, if not all, should 
be resolvable. Besides TON 139 and PG 1718+519 we 
find in three cases companions to the sdB stars which are 
too faint to match the spectral energy distribution. These 
could be triple systems consisting of an unresolved sdB 
binary and a distant third star. 

6. Radial velocities 

Important additional information can be obtained from 
radial velocity measurements. A systematic search for ra- 
dial velocity variations of our programme stars is needed. 
Such projects have alread y bee n started by Saffer et al. 
(2001) and Maxted et al. (2001) who observed six of our 
programme stars (PB 6107, PHL 1079, PG 0749+658, 
TON 1281, PG 1449+653 and PG 2148+095). None of 
them showed significant radial velocity changes. 

Saffer et al. ( ^00l| ) find in their survey of 21 compos- 
ite spectrum sdB stars that the velocity variations of the 
individual components as well as the velocity difference 
between the two components are very small (less than a 
few kms^^) or undetectable, and conclude that the bina- 
ries have likely periods of many months to several years. 
Green et al. (2001) estimate from these measurements that 
the current periods average 3-4 years with separations 
540 - 650 Rq. 

We have obtained multiple precise radial velocities for 
TON 139 and a single measurement of PG 1718+519 us- 
ing the MMT Blue Channel spectrograph at lA resolu- 
tion from 4000-4930 A (see table The radial velocities 
of the cool companions were determined by cross corre- 
lation against super-templates of main sequence spectral 
types from F6 to K5. The sdB velocities were derived us- 
ing a preliminary attempt at subtracting out the cool star 
companion spectrum. For details of the data reduction 



^ Note that PG 0105+276, which is resolved in three com- 
ponents (see Fig. is an sdO star and does not belong to our 
sdB sample. 



and analysis see Saffer, Green & Bowers (2001). Improved 
sdB velocities using better cool star template spectra for 
the subtractions will be determined by Green, Bowers & 
Saffer (2002, in prep.). 

For TON 139 the cool star's velocity is constant, 
whereas the sdB velocity is changing by more than 
50kms^^. This can be explained if an additional com- 
panion is orbiting the sdB star. This companion has to be 
so faint that it does not contribute to the light in the R 
band. Hence we have to conclude that the resolved sys- 
tem TON 139 is a triple system. A radial velocity study 
of PG 1718+519, the second resolved system in our sam- 
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Table 6. Heliocentric radial velocities for the sdB- and the cool star components of TON 139 and PG 1718+519. 
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pie, is not available yet. The single measurement listed in 
Table ^ gives identical radial velocities for the sdB and the 
cool companion. This argues against a third faint compo- 
nent orbiting the sdB star in a narrow orbit as was found 
for TON 139. Additional radial velocity measurements are 
urgently needed to clarify the nature of PG 1718-f519. 
Assuming that PG 1718-1-519 is not triple, this would be 
the only resolved binary system in our sample of 19 ob- 
jects. 

7. Conclusions 

In total we have resolved six systems out of a sample of 23 
stars. Of those 23 stars, however, four do not really belong 
to the intended sample of sdB stars showing evidence for 
a cool companion: PG 1558-007, KPD 2215+5037, and 
PG 2259+134 show no photometric or spectroscopic evi- 
dence for a companion. The observed infrared excess can 
be explained by interstellar reddening rather than by a 
cool companion. 

PG 1558—007 does have a resolved near by star (linear 
separation 1500 AU), which, however, is too faint to con- 
tribute detectably to the combined light in the R band. 
PG 0105+276 is a helium-rich sdO star (with two possible 
distant companions at 3700 AU and 4900 AU). 

Of the remaining four resolved systems the nearby 
stars are in two case (TON 1281, HE 0430-2457) too faint 
to reproduce the photometric and/or spectroscopic obser- 
vations of the stars. 

Only in the two systems TON 139 and PG 1718+519 
(separations 0'.'32 and 0'.'24, respectively) do the magni- 
tudes of the resolved components match the expectations. 
These two stars could be physical binaries whereas in the 
other cases the nearby star may be a third component or a 
chance projection. Radial velocity measurements indicate, 
however, that the resolved system TON 139 is also triple. 

Hence, the observed sdB binary sample was reduced 
to 19 objects with two bona-fide resolved systems, which 
have apparent separations of 0'.'24 and 0'.'32. From the nu- 
merical simulations we would expect to resolve six to seven 
systems if sdB stars have the same binary characteristics 
as normal stars, out of which one system is expected to 
have a sin I >1" and two should have separations between 
O'.'l and 0'.'2. The discrepancy becomes even more pro- 
nounced if one recalls that our photometric fit procedure 



tends to underestimate the brightness of the companion 
(and thus to overestimate the limiting angular separation 
that can still be resolved). In addition we expect three 
triple systems to be present in our sample. Most of these, 
if not all, should be resolvable. Such systems could explain 
some of the more distant companions as well as the radial 
velocity measurements of TON 139. 

This success rate (1 resolved binary out of 19 candi- 
dates) is clearly below the prediction of numerical simu- 
lations assuming single star evolution (about 30%), using 
the distribution of binary separations given by Duquennoy 
& Mayor (1991). This indicates that the distribution of 
separations of sdB binaries strongly deviates from that of 
normal stars. 

If, on the other hand, all sdB stars were produced by 
close binary evolution, none of the binary systems should 
have been resolved (even at the high spatial resolution of 
the WFPC2 camera). Our low success rate is thus closer to 
that predicted by the close binary evolutionary scenario. 



Recent radial velocity surveys (Saffer et al. 2001; Maxted 
et al. 2001) revealed that a large fraction of single- lined 
sdB stars are indeed close binaries with periods below 10 
days. Our results could be explained if most of the pro- 
gramme stars were close binaries. Therefore, our study 
provides further evidence that close binary evolution in- 
deed is fundamental to the evolution of sdB stars. A survey 
for radial velocity variations in all of our programme stars 
will be tale telling. 
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Appendix A: Spectroscopic observations and data 
reduction 

The observational setups and observing dates for the new 
spectra are given in Table A.l. The reduction of the 
spectra of PG 0105+276, HE 0430-2457, PG 0942+461, 
HE 2213-2212, and K PD 2215+5037 are described by 
Edelmann et al. ( ^001b| ). PG 2148+095 was o bserved and 
reduced as described by de Boer et al. ( 1995| ), the reduc- 
tion of PG 1309-078 and PG 1558-007 was perfor med in 
the same way as described in Moehler et al. (1997). 

Fig. [A.2| shows the spectra of the stars that show no 
spectroscopic or photometric evidence for a cool compan- 
ion (PG 1558-087, KPD 2215+5037, and PG 0105+276). 
The Ca ii absorption lines in the spectra of these stars (see 



Fig. A. 2) are probably of interstellar nature. Our spectrum 
clearly shows that PG 0105+276 is a helium rich sdO star 



(see Fig. A. 2) inconsistent with the photometric classifi- 
cation as sdB+K7 by AUard et al. (1994, where all three 
stars seen in Fig. ^ were included in the measurements) 
but in accordance with the early spectroscopic classifica- 
tion by Green et al. ( |1986| ). 

We derived the atmospheric parameters Toff, log g and 
helium abundance simultaneously for the single stars by 
matching a grid of synthetic spectra derived from H and 
He li ne blanketed NLTE model atmospheres (Napiwotzki 



1997) to the data. For temperatures below 27,000 K we 



used the met al lin e blanketed LTE model atmospheres of 
Heber et al. (200C). The synthetic spectra were convolved 
beforehand with a Gaussian profile of the appropriate 
FWHM to account for the instrumental profile. Results 
are given in Table A.l and Fig. A. 3 displays the fit for 
KPD 2215+5037 as an example. 




x/k 



Fig. A.l. Comparison of normalized spectra of four pro- 
gramme stars to PG 1309—078, which is already known to 
be a spectroscopic binary containing an sdB. The spectral 
features indicative of a cool companion are marked. 





He I He I 




Fig. A. 2. Optical spectra of the sdB stars PG 1558-007 
and KPD 2215+5037 as weU as of the sdO star 
PG 0105+276. The spectra of the former are dominated 
by hydrogen lines whereas that of the latter by He ii lines. 




Fig. A.3. Spectral fit for the sdB star KPD 2215+5037. 
He is excluded from the fit because of contamination by 
interstellar Can. 

Appendix B: Photometric data for our 
programme stars 



In Tables B.l and B.2 we compile the photometric data 



collected from literature and used in the photometric de- 
convolution. 
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Table A.l. New optical spectroscopy and atmospheric parameters of single programme stars 
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star 


telescope and 


wavelength 


spectral 


obs. date 




log g 


log(He/H) 




spectrograph 


range resolution 














[A] 


[A] 




[K] 


[cgs] 




PG 0105+276 


CA 3.5m TWIN 


3600 - 7400 


3.1 


1997/08/31 


63000 


5.4 


+0.5 


HE 0430-2457 


ESQ 1.5m B&C 


3600 - 7450 


5.5 


1996/10/22 








PG 0942+461 


CA 3.5m B&C 


3860 - 5560 


5.0 


1989/01/23 








PG 1309-078 


ESQ 1.5m DFOSC 


3860 - 6780 


5.4 


2000/06/21 








PG 1558-007 


ESQ 1.5m DFOSC 


3860 - 6780 


5.4 


2000/06/21 


20300 


5.0 


-2.6 


PG 2148+095 


ESQ 1.5m B&C 


3730 - 4970 


3.0 


1991/07/10-15 








HE 2213-2212 


ESQ 1.5m B&C 


3600 - 7400 


5.5 


1996/10/23 








KPD 2215+5037 


CA 3.5m TWIN 


3260 - 7450 


3.1 


1997/08/29 


29400 


5.6 


-2.2 


PG 2259+134 




Theissen et al. 


1993 




31900 


5.9 


-1.7 



Table B.l. Stromgren photometry and UV spectrophotometry for our programme stars. Stromgren photometry is 



taken from Gree n (pM3| , G80), Kilkenny ( p^84[ K84;|l98^, K87), Moehlcr et al. ( |l990| , M 90), Theissen et al. ([19931, T93 ), 
Wesemael et al. (1992, W92). The lUE data were obtained from the lUE final archive (http://archive.stsci.edu/iue/). 



Star 


y 


b-y 


u — b 


nil 


Cl 




Ref. 


lUE 

SWP 


LWP 


PB 6107 


12'^897 


+0"'032 


+0^112 


+0'"052 






W92 








12'^889 


+0'"026 




+0"092 


-0 


"094 


M90 








12"89 


+0'"01 


+0'^'10 


+0™05 






G80 








12™907 


+0'"018 




+0"093 


-0 


"109 


K87 






PG 0105+276 


14™481 


+0™022 


-0^194 


+0'"023 






W92 


56271 




PHL 1079 


13™278 


+0"^'003 


+0"^106 




-0 


■^109 


K84 


42338 


21098 


PG 0749+658 


12^135 


-0^032 


+0^131 


+0"^087 






W92 






TON 1281 


13'^371 


+0"'094 


+0^175 


+0'^'065 






W92 


56384 




TON 139 


12™796 


+0™111 


+0^^364 


+0"055 






W92 






PG 1309-078 


14'"11 


+0'"07 


+0'^'06 


+0"18 






G80 






PG 1449+653 


13^580 


+0™041 


+0'^'047 


+0"034 






W92 


34298 




PG 1511+624 


14™421 


+0";'049 


-0^002 


+0'^'005 






W92 


39370, 57359, 57361 


18491 


PG 1558-007 


13'^528 


-0^011 


+0™244 


+0'"091 






W92 






PG 1636+104 


14'"090 


+0"'169 


+0™426 


+0'^056 






W92 






PG 1656+213 
















39422 


18542 


TON 264 


14'^'070 


+0"'008 


-0"^053 


+0'^070 






W92 


39422 


18542 


PG 1718+519 


13^686 


+0^102 


+0"^307 


+0"084 






W92 


41571 


20308 




13*^694 


+0™131 




+0^094 


-0 


■^095 


T93 






PG 2148+095 


13'^037 


+0"'028 


+0"^087 


+0'"066 






W92 


56148 




KPD 2215+5037 


13'^739 


-0™026 


+0^034 


+0'^068 






W92 






PG 2259+134 


14™478 


-0'"038 




+0^082 


-0 


"^089 


M90 


44821, 56182 


23244 


PG 2259+134 


14"^545 


-0'"069 


-0™011 


+0^088 






W92 






BD -7°5977 
















31030 


10815 
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Table B.2. BVRI (AUard et al. [l994| ), 
frared broadband photometry (UT98: UUa 



UBVI (F ergus on et 
& 



ThejU 



1998 



http://irsa.ipac.caltech.edu/applications/2MASS/BasicSearch/) for our programme stars 



al. 

2MASS: 



1984| ), HST R (this paper) and in- 
2MASS 2"^^ incremental data release, 





V 


1 ) — V 


T/ R 
V — iX 


It — 1 


-Rhst 


J 

J 


±1 


1^ 
j\ 


Pof 
IXCi. 


rSD DlU/ 


Iz • ooi 


— U • Uoo 


+U • U /U 


+u ■ uyo 


Iz • oU 










■pr^ ninp;_i_07'fi 
r\j UlUO+Z/D 


14 ■ 44o 


— U • Uo / 


+U ■ UoD 


+U ■ Iz / 


14 • OD 


1 A 7 
14 • o4 / 


lo • ozl 


lo • / zl 


OA/T A GQ 

zMAdo 


T")T TT -1 r\'-7r\ 










lo • z4 


Iz • Oo 


Iz • ZO 


Iz • U4 


u 1 yo 


TTTT ri/IQA 0/1 
rirj U4oU — Z40 ( 


14 • 100 


— U • U4D 


+U • UoO 




1 A "^07 

14 • U ( 


lo • Diy 


lo • olO 


lo • zUo 


OA/T A QQ 

ziVlAoo 




1 o™i oi 
Iz ■ Izl 


o"^i O^i 
— U • lUD 


+U ■ Uzl 


+U ■ U / z 


Iz • 14 










1")/' 1 nri/io 1 /ifii 










lo • yo 


lo • DlZ 


lo • 1 / Z 


lo • Uo4 


OA/T A QQ 


TOM lOfil 


lo ■ 4oy 


-f-u • Uy4 


-f-U ■ lOD 


1 n"^i 7fi 


1 q"^07 
lo • z / 


Iz ■ / Oo 


Iz • OUo 


1 o"^A A S 
Iz • 44o 


9A/T A QQ 


IVJIN ioy 










Iz • DO 


1 o™i n 
iz • lu 


11 • yz 


11 • yo 


u 1 yo 


r^Vjr loUy — U/O 










14 • UO 


lo • OOo 


lo • zoy 


lo • IDz 


OA/T A QQ 


PG 1449+653 


13^^611 


-0^^035 


+0'"073 


+0^110 


13™57 










1 V.I J_<J J_ T-p-UZ^ 


±^ OZ ( 


— u uzz 


1 n™i 1 


1 n™i AO 


1 zt™^8 




TO OlO 


±0 OOO 


OA/T A 


1 V.I lOOo — UU 1 


lO ■ 041 


0"^0fi/1 
— U • U04 


^n"^m 9 

• UiZ 


-f-U ■ IIU 


1 Q^^Kc: 
lo • 00 










pr^ 1 fioi _Li /I 


14 ■ 4z4 


_l0"^092 
-f-U • Uzo 


-f-u ■ ioU 


-f-U ■ o4 ( 


1 /I "^Q7 


10 • yi(5 


lo • i o 


lo • DUU 


OA/T A QQ 
ZiVlAoo 


Pr^ 1fiQfi_LlO/1 


1 A ™nQQ 

14 ■ uoy 


-)-u • lyo 


-f-u • iyi 


+U • lyD 


1 q'^qp: 
lo • oO 










TON 264 


14^074 


—0^^083 


+0™066 


+0™136 


14™02 










PG 1718+519 


13^^733 


+0'^113 


+0"'156 


+0"132 


13"53 


13^^008 


12^^716 


12'^664 


2MASS 


PG 2148+095 


13^^021 


-0'^024 


+0"'060 


+0'"096 


12'"98 


12"'18 


12^^34 


12'^06 


UT98 


HE 2213-2212 










M-^Ol 


13"'686 


13^^292 


13'^236 


2MASS 


KPD 2215+5037 


13'^664 


-0^^093 


+0"'015 


+0'^052 


13^91 










PG 2259+134 










14^70 


15"'795 


15^^220 


14^^523 


2MASS 


BD -7°5977 


10^55^ 


+0^51' 






10'"05 


8'^97 
9^^017 


8^M8 
8"'526 


8'^38 
8^^448 


UT98 
2MASS 


Star 


V 


B-V 


U-B 


V -I 


-Rhst 


J 


H 


K 


Ref. 


PG 1421+345 


14^^78 


-0^14 


-0'"89 


+0'^63 


14^59 


14"'035 


13^^716 


13'^678 


2MASS 


PG 1601+145 


14^^50 


+0'^01 


-0™92 


+0'"41 












PG 1656+213 


14^88 


-0'^20 


-0'"73 




14"'73 











^ Altmann (priv. comm.) 

^ Derived from Tycho photometry (Vr 



10 -6, (B — V)t = +0 • 6) using the transformation given in Ferryman (1997) 



